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Abstract 

We study the steady state structure and dynamics of a 2-d Ising interface placed in 
an inhomogeneous external field with a sigmoidal profile which moves with velocity 
Ve- In the strong coupling limit the problem maps onto an assymmetric exclusion 
process involving motion of particles in 1-d with position dependent right and left 
jump probabilities. For small Ve, the interface is stuck to the field profile. As 
increases the profile detaches from the interface. At the transition point( and beyond 
), the interfacial structure and dynamics is characterized by KPZ exponents. For 
small Ve, on the other hand, the interface is macroscopically smooth with a vanishing 
roughness exponent a. The interfacial structure is periodic with a periodicity which 
depends on the orientation of the interface. For a fixed orientation this periodic 
structure "melts" as Ve is increased. We determine the dynamical "phase - diagram" 
of this system in the Ve - orientation plane. 
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1 Introduction 



Consider a one- dimensional (1-d) Ising interface [1,2] between up and down 
spins in two dimensions (2-d) obeying single-spin flip Glauber dynamics [3,4]. 
In the presence of a uniform driving field [5-10] the interface moves with a 
velocity which depends on the magnitude (and sign) of the driving field. On 
the other hand, a fixed external field profile which is positive in the region 
of up spins and negative in the region of down spins, would stabilize a sta- 
tionary, macroscopically fiat interface. In this paper we study systematically 
the structure and dynamics of this Ising interface as this field profile is moved 
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Fig. 1. An Ising interface u{y,t) (bold curved line) between regions of positive 

(marked +) and negative (marked — ) magnetization in an external, inhomogeneous 
field with a profile which is as shown(dashed line). The positions of the edge of the 
field profile and that of the front are labelled Se and Sf respectively. 

with an arbitrary velocity Ve- 

There are two reasons why we are interested in this problem : Firstly, there 
are several important practical examples where inhomogeneous fields drive 
interfaces. Some of them include zone purification of Si where the controlled 
motion of a temperature field profile is used to preferentially segregate im- 
purities[ll], magnetization of a bar of iron with a permanent magnet, phase 
transitions induced by a travelling heat (welding) or pressure (metamorphosis 
of rocks) fronts etc. Secondly, we would like to extend this study, in the fu- 
ture, to the dynamics of solid interfaces where interfacial degrees of freedom 
are coupled to hydrodynamic modes of the bulk solid eg. phonons (acoustic 
emmissions[ll]) and defects[12]. A systematic study of how these modes are 
excited in sequence as is increased is of great fundamental interest. 
The object of our study here is the interface between up and down spin phases 
(Fig. 1) in the limit h/J,T/J —>■ 0, where J is the Ising exchange coupling, 
T the temperature and the field h{x, t) here is inhomogeneous, h — hmax in 
the region where the magnetization is positive and —h^ax in region where it is 
negative separated by a relatively sharp edge. The edge of the field (i.e. where 
the field changes sign) lies at S^. The front or interface, u(y,t) (no overhangs 
!) separates up and down spin phases. The interface is a bold curved line with 
the average position Sf. To move the interface we move the edge with velocity 
Ve] in response the front moves with velocity Vf. Parts of the front which leads 
(lags) the edge of the field experience a backward (forward) force pulling it 
towards the edge. The driving force therefore varies in both space and time 
and depends on the relative position of the front compared to that of the edge 
of the dragging field. What is the behaviour of the front velocity Vf as a, func- 
tion of Ve? What is the structure of the interface in various regimes? These 
are the questions we address in this paper. 

Briefiy, our results are as follows 

(1) For any orientation of the interface, Vf = for small the front moves 
along with the field profile. We call this the "stuck" phase. 

(2) At a velocity Vg = v^, the front detaches from the field profile. At higher 
values of Ve the front experiences an uniform magnetic field h — h^ax 
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and the problem reduces to the growth of an Ising interface driven by a 
uniform field [5-10] 
(3) The structure of the "stuck" interface is flat with a roughness exponent 
a — 0. Depending on the orientation of the interface, the height of the 
interface u{y,t), as a function of y and time t, may show periodic oscil- 
lations in y and/or t. The nature of these oscillations depends crucially 
on the system size in a manner to be explained below. 

In the next section we introduce a continuum description of the problem and 
derive the relevant coarse grained equations of motion. In Section 3 we present 
the mean field solution to these equations. In Section 4 we introduce fluctu- 
ations through an exact mapping to an assymmetric exclusion (particle hop- 
ping) model in 1-d and present, analyze and discuss our results obtained from 
computer simulations of this model. In Section 5 we present our conclusions. 



2 Continuum Description 

Let the magnetization of the 2-d Ising system be given by, (p = (j){x,y,t). We 
assume that for h/ J,T/ J ^ the magnetization is uniform everywhere except 
near the interface which may be parametrized by a function u{y,t), where u 
is the height of the interface perpendicular to y. Hence the magnetization 
= (f){x — u{y,t)). The field profile is given hy h = hmaxis.nh.{{x — Vet)/x) 
where x is the width of the profile. Model A dynamics[13] for then implies, 

where 

HT = j dr[ai<j)^ + a2<f)^ + asiWcpf - h{x, t)<f)\ (2) 

is the Hamiltonian and C is a Gaussian white noise with zero mean and 

< C(r, t)C{r\ t') > = 2kBTr5{r - r')5{t - t') (3) 

Using Ht in Eq.l, taking x — u{y,t) = v and converting all derivatives to 
derivatives over the profile u{y,t), we have. 
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- ct>\vf-^ = - T[2a,ct>{v) + Aa2ct>\v) - 2a^4>\v) + 2a30'(t;)^ 

-2a,ct>"{v){-f-h{x,t)\-rCM (4) 



We then choose a (f) dependent mobihty F contributing to the lowest order 
in (f) consistent with symmetry viz. F = Fq + Fi(V0)^. Substituting for F 
and integrating both sides of the equation with respect to x between hmits 
{u — x/2) and {u + x/2) i-e. over the interfacial region, remembering that 
has a sigmoidal profile, we finally get an equation of motion for the profile u. 



du , ff^U , /du\2 /U — Vet\ , , /U — Vet\ , 

- = A,— -A,(-) tanh(— -^)-A3tanh(— -^)+C'k^) (5) 



where Ai,A2 and A3 are constants. This is different from the familiar KPZ 
equation [10] in the fact that it lacks Galilean invariance { u' ^ u + ey, — > 
y - X2et, t' ^t) 

In general crystal field effects introduce a lattice periodic force [1] which may 
be accounted for by including an additional term Vq sm{27Tu/a) (a is the lattice 
parameter) to the right side of the above equation of motion. 



3 Mean Field Result 



A mean field calculaton amounts to taking u = u{t) i.e. neglecting spatial 
fiuctuations of the interface. Then 



du , , /U — Vet\ . 

— = - A3 tanh (-^) + C («, t) (6) 



For large times {t 00), u Vft, where Vf is the average velocity of the 
front. Thus Vf is obtained by solving the self-consistency equation, 

., = -A3tanh(i^^^^) (7) 



In the t ^ 00 limit the hyperbolic tangent is replaced by the simply the sign of 
Vf — Vf. and is equivalent to taking x — > namely, an infinitely sharp field pro- 
file. For small Ve the only solution to the self-consistency equation is = I'e 
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as can easily be verified graphically. For large edge velocities > v^, where 
V* = A3 wc get f J = A3 = V*. We thus have a sharp transition from a region 
where the front follows the edge with the same velocity to one where it moves 
with a constant velocity unable to follow it anymore. The interface velocity 
relaxes to its steady state value Ve as 1/t in the region of low Ve- The region 
where the front moves with a constant velocity is evidently the well studied 
problem of driving an interface by a homogeneous field [6,9]. 
How is this result altered by including spatial fluctuations of m ? In order 
to answer this question we have mapped this interface model to an assym- 
metric exclusion process [14] and study the dynamics both analytically and 
numerically using computer simulations. 



4 Beyond Mean Field Theory 

The mapping to the exclusion process follows [8] by considering Np parti- 
cles distributed among Ng sites of a 1-d lattice. The particles are labelled 

n = 1,2, ,Np sequentially at t = 0. Any configuration of the system is 

specified by the set of integers {y{n)} where y{n) denotes the location of the 
nth particle. In the interface picture n needs to be interpreted as a horizontal 
coordinate (y in Fig. 1), and y{ri) as a local height u{y,t). Each configura- 
tion {y{n)} then defines a one-dimensional interface inclined to the horizontal 
with mean slope 1/p where p = Np/Ns. The interface coordinates satisfy 
y{n + 1) > y{n) + 1, and periodic boundary conditions amount to setting 
y{n + Np) = y{n) ± Ng. Motion of the interface under the influence of a 
driving field corresponds to the hopping of particles. In each time step {Np 
attempted hops), y{n) tends to increase (or decrease) by 1 with probability 
p (or q); it actually increses (or decreases) if and only if |/(n + 1) — y{n) > 1. 
In our case the right and left jump probabilities p and q {p + q = 1) are 
not constants but themselves depend on the relative position of the inter- 
face y{n) and the edge of the field n/p + Vet. Note that in calculating this 
relative position we have to use the actual position of the interface without 
periodic boundary conditions. In our calculations reported here we use a bias 
A = p — q = Aosign(|/(n) — n/p — Vet) with Aq = 1 unless otherwise stated. 
We are interested in the average vertical velocity of the interface Vf defined as 
the total number of particles moving right per time step. In addition to the 
front velocity, we also examine the behaviour of the width of the interface: 

<y\t) = Np' E <{y{n,t)-y{n,Q)-Vftf> (8) 

n=l,Np 

as a function of time and system size Ng. The angular brackets denotes an 
average over the realizations of the random noise. 
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Fig. 2. The front velocity f/ as a function of the velocity of the dragging edge 
for = lOO(n), 1000(0), 10000(+) and p = 0.5. All the data (symbols,n, O, +) 
collapse on the mean field solution (dashed line). 

Note that the usual particle hole symmetry for an exclusion process [5] is 
violated since exchanging particle and holes changes the relative position of 
the interface compared to the edge of the field. 

J,..! Monte Carlo simulations for the dynamical transition 

We perform Monte Carlo simulations of the exclusion process using a ran- 
dom sequential update [15] to understand the behavior of the interface in the 
presence of the inhomogeneous field profile. We study the system for different 
system sizes and densities and obtain the average velocity of the steady 
state interface as a function of the velocity of the edge of the field profile v^.. 
The velocity Vf is obtained by dividing the distance moved by the interface for 
a certain (large) number of time steps by the total number of time steps after 
discarding the first few thousand steps to remove transients. Fig. 2 shows a 
sharp dynamical transition from an initially "stuck" interface with Vf = Veio 
a free, detached interface with Vf — vl — — p) the result for an assymet- 
ric exclusion process with density p. Note that the mean field solution for the 
front velocity and the dynamical transition is exact. 

4-2 The Stuck phase (vg < v^) 

The stuck phase is characterized hy Vf — Vg and a bounded. To obtain the 
ground state of the interface in the presence of a stationary = 0) field 
profile it is sufficient to minimize I]„(y(n) — n/ pY which demands y{n) to 
follow the edge n/ p as closely as possible subject to the constraint that y{n) 
be an integer. This ground state structure is always periodic for p < 1/2. For 
p > 1/2 this periodicity is destroyed for infinetisimal v^- For densities which 
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are of the form \/k where k is an integer, the result is particularly simple 
viz. y{n) = kn. This corresponds to a particle-hole system where the particles 
form a 1-d lattice with a lattice parameter of k. For an arbitrary density 
(orientation) the ground state is still periodic over short distances but has long- 
period (possibly incommensurate) modulations. We verify this by calculating 
the pair distribution function g{l) — {l/Np{Np — l))J2n,m^i,{y{m)-y(n)), where 
n,m are particle indices and m > n. The Fourier transform of g{l) shows 
prominent delta function peaks. 

The dynamics of the interface for Ve > depends on whether or not the system 
size is compatible with the lattice parameter k. If the system size A^^ is an exact 
multiple of k then the particles which are separated by intervening regions of 
holes can move independently of each other in response to the local bias A. 
Let Pi be the probability {J2iZ-oo — 1) obtaining a particle (any particle) 
in state i where i — Q corresponds to a particle which has not moved from its 
initial position and i — s {— —s) corresponds to one which has moved s integral 
lattice spacings to the right (left) of its original position. Fluctuations of the 
interface about the ground state correspond to these independent particle 
motions which cost energy if v^t is integral. Consider now that i < v^t < i + \ 
the form of the bias A implies that the probabilities satisfy the following 
set of master equations, 



P. 



-Pj-i ~ Pj + P^ 



-P, + P. 



for j > i + 1 
for j = i, i -|- 1 
for j < i. 



(9) 



Noting that the average position S{t) =< N^ ^ Uni.uiP'ii) ~ uij'')^) > the 
interface is given simply by S{t) = J2iZ-oo Pii't) and a^{t) = Ei^-oo(^^ ~ 
i) Pi{t) we obtain the results shown in Fig. 3(a) and 3(b). It is clear that these 
results match the corresponding ones obtained from Monte Carlo simulations 
exactly. The interface therefore follows the profile in a jerky fashion and the 
width of this interface oscillates between fixed bounds. Thus although the 
structure of the moving interface corresponds more or less with the ground 
state periodic structure for small v^i the velocity is oscillatory. If however p is 
not the reciprocal of an integer and the system size does not accomodate an 
integral niimbcr of spatial periods of the ground state then the corresponding 
1-d lattice contains long wavelength modulations and the particle hoppings 
are not independent anymore. For fixed edge velocities it is found that is a 
constant in time independent of the size of the system, the constant, however, 
depends on and p. The average position of the interface does not show any 
oscillations and faithfully follows the field profile corresponding exactly to the 
mean-field solution. 

The Melting Transition 
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Fig. 3. (a)Variation of with t for for p = 0.5, Vf> = 0.025 and p = 1.0 and (b) 
variation of S{t) with t for p = 0.5, = 0.025 and p = 1.0. Lines denote analytical 
results while points denote monte carlo data. 

The Fourier transform of i.e. the structure factor g{q), indicates a "melt- 
ing transition" of the periodic steady states with increasing Ve for any density 
(p < 1/2). We track this by plotting the intensity of the largest peak (smallest 
q) of g{qmax) as a function of Vg for a number of densities. This is shown in 
Fig. 4(a) . Also, the lattice paramater (a) obtained from 27r /q^ax is used to 
determine the Lindemann ratio L = o"^ jd^ . The increase of the Lindemann 
ratio with the fg for different densities (Fig. 4(b)) is another proof of a melt- 
ing transition. These results are summarized in the dynamical phase diagram 
(Fig. 5) for the Ising interface in a moving field profile. It shows an detachment 
transition along the line = (1 — p) (for Aq = 1) and a melting transition. 
The exact position of this dynamical melting transition (unlike a thermody- 
namic transition) depends on the parameter used to characterize it. If we use 
giflmax) then the melting transition occurs simultaneously with detachment 
for p < 0.5 and at Vf. = for larger p. Using the Lindemann parameter, 
however, one obtains a melting transition which preempts detachment. 



5 Behavior at the transition point 

We want to determine scaling form for a{t) at the transition point viz. the 
growth exponent /?, the roughness exponent a and the dynamic exponent z. 
In the detached phase we know from renormalization group analysis that the 
exponents are in the Kardar-Parisi-Zhang (KPZ) universality class [1,10] viz. 
P — 1/3, a = 1/2 and z — a/P — 3/2. To determine these exponents at the 
transition point we make use of Family-Vicsek scaling relation [1] 

aiL,t)r.N:f{t/N:) (10) 

Fig. 6 shows the variation oit/N^ with a{L, t) /Nf for different p and different 
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Fig. 4. (a) The structure factor g{qmax) vs. v^, for various p. Note that g{q,nax) 
vanishes as Vg increases thereby implying a melting transition in the 1-d assymetric 
exclusion process, (b) Variation of lindemann ratio L with edge velocity Ve for 
various densities. An increase of L again signifies melting. 
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Fig. 5. Phase diagram for Ising interface driven by an inhomogeneous magnetic field 

showing periodic stuck (PS) steady states for < (1 — p) and p < .5, uniform stuck 
(US) steady states for Ve < {I — p) and p > .5 and detached (KPZ) steady states 
for fe > (1 — p)- 

system sizes A^^. The curves collapse onto one curve once an intrinsic width 
CTj, arising from finite-size and crossover effects [1], is substracted out. The 
exponents were found to be KPZ. 

To understand why this happens we go back to our modified KPZ equation 
(Eq.5) and make the transformation u = u' + Vft. We get 



9i/ 
'dt 



+ ^;/ = Al— -- Asl^ tanh(^-^ 

{Vf — Ve)t + u'-' 



Ve)t + m'n 



—A3 tanh ^- 



X 



-) + C{u,t) 



(11) 
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Fig. 6. Monte Carlo data {a - ai)/Nf vs t/N^ for p = 1.0, p = 0.7 and with 
Ns = 1000 [x,0] , 800 [0,0] , 400 [+,*] , 200 [O, A]. All the curves collapse to a 
single universal function showing KPZ scaling. 

Now substituting the mean field result for Vf (Eq.7), making use of the fact 
that at the transition point = A3 and simplifying one can show that the 
above equation reduces to the famihar KPZ equation in u' . 



6 Conclusion 



In this paper we have studied the static and dynamical properties of an Ising 
interface in 2-d subject to a non-uniform, time-dependent external magnetic 
field. The system has a rich dynamical phase diagram with several dynami- 
cal phases (steady states). The nature of these steady states depend on the 
orientation of the interface and the velocity of the external field profile. The 
detailed dynamics of the interface depends on whether or not the size of the 
system is commensurate with the orientation. For a commensurate system, 
the interface follows the field in a jerky fashion and the width of the inter- 
face fluctuates between well defined bounds. For a general, incommensurate 
interface the motion of the interface is steady and the width is constant. For 
large velocities of the external field, the interface detaches from the profile and 
coarsens over time with KPZ exponents 

In future we would like to study in detail further dynamical aspects of this 
system e.g. the hysteretic response of this system under time varying external 
parameters {ve). 
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